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1 .  RESEARCH  OBJECTIVES 


(a)  Combustion  Instability 


Over  the  past  five  years  we  have  studied  the  detailed  mechanism 
of  combustion  instability  in  a  simple  dump  combustor.  Our 
determination  to  understand  the  mechanism,  rather  than  the  symptoms, 
was  driven  not  only  by  the  -tact  that  other  such  studies  seldom  it 
ever  get  to  the  mechanism  but,  toward  the  completion  of  the  work,  by 
the  requirement  to  show  that  efficient  and  simple  control  of 
non-linear  instability  reouires  a  knowledge  of  the  mechanism. 

The  effort  this  year  was  devoted  i)  complete  studies  of 
chemical  time  delay  and  ii>  to  carry  out  spatial  heat 
re  1  ease-pressure  correlations  showing  the  manner  in  which  the 
processes  within  the  entire  chamber  developed  into  the  limit  cycle. 
The  first  constituted  the  completion  of  work,  a  portion  of  which  was 
reported  last  year.  The  second  was  a  new  and  extensive  set  of 
measurements  which  were  of  particular  interest  in  demonstrating  the 
essential  non-linearity  of  the  mechanism. 

This  will  complete  our  work  on  the  mechanism  of  combustion 
instability;  any  further  work  will  deal  with  the  control  of 
combustion  instability  and  will  be  the  subject  of  a  new  proposal. 


< b)  Experiments  on  Combustion  in  Large  Vortices 

The  primary  thrust  of  this  grant  concerns  a  detailed  study  into 
the  combustion  process  of  large  vortices.  It  is  closely  related  to 
our  extensive  effort  on  the  possibilities  of  shock  enhancement  of 
mixing  and  combustion  of  hydrogen  with  air. 

The  principal  effort  this  year  has  been  to  design  and  build  a 
new  combustion  chamber  which  will  allow  the  production  of  large 
vortices  in  a  we  1 1 -con tr ol 1 ed  periodic  manner  and  over  a  wide  and 
controllable  range  of  frequencies.  At  the  same  time,  we  are  making 
major  improvements  in  the  data  acquisition  system,  centering  about  a 
CID  array  camera  with  image  intensification  and  an  LDV  for  the 
purpose  of  making  detailed  velocity  measurements  in  the  hot  gas. 
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<  c )  Theoretical  Studies  of  Cortex  Combustion 

The  theoretical  investigation  ot  flames  distorted  by  vortex 
structures  was  initiated  under  this  grant  some  years  ago  and  has 
proven  to  be  an  extremely  informative,  as  well  as  popular,  direction 
of  research.  With  regard  to  the  present  program,  the  aim  of  these 
investigations  is  both  to  provide  i n ter pre tat i on  of  the  experiments 
and  to  develop  laws  that  allow  scaling  the  results  to  different 
parametric  ranges.  The  work  this  year  yielded  important  information 
concerning  the  mixing  and  chemical  reaction  when  the  viscous  effects 
were  1 arqe  enough  to  retard  the  deformation  of  the  flame  and  hence 
to  decrease  the  overall  combustion  rate.  A  new  dimensional  quantity, 
the  ratio  flame  sheet  spreading  to  viscous  diffusion,  was  introduced 
as  a  measure  of  this  phenomenon. 


2.  STATUS  OF  RESEARCH 


<  a)  Combustion  Instability 


An  investigation  of  the  mechanism  leading  to  growth  and 
sustenance  o-f  large  pressure  oscillations  in  a  laboratory  combustor 
has  been  completed.  This  work  has  shown  that  relatively  low 
■frequency  oscillations  can  be  produced  in  airbreathing  systems  o-f 
reasonable  size  and  that  these  oscillations  are  associated  with 
longitudinal  acoustic  modes  of  the  entire  combustion  duct.  The 
features  of  this  mechanism  have  been  made  clear  by  this 
investigation  and  it  is  sufficiently  understood  that  we  can  now 
develop  methods  for  its  suppression. 

The  mechanism  which  underlies  the  development  and  sustenance 
of  the  pressure  oscillations  contains  the  following  elements.  Ule 
consider  a  combustor  supplied  with  a  fuel-air  mixture  in  which  the 
flame  is  stabilized  on  a  rearward  facing  step  as  shown  in  Fig.  1. 
The  mechanism  is  the  same  for  systems  in  which  bluff  body  flame 
holders  such  as  “U“  gutters  are  used. 


1.  Combustion  Oscillations  on  Rearward  Facing  Step. 

Pressure  oscillations  at  one  of  the  longitudinal  acoustic  modes 
of  the  system  produce  velocity  and  pressure  fluctuations  at  the 
flame  holder  lip  which  periodically  trigger  the  formation  of  large 
vortices  in  the  recirculation  zone  downstream  of  the  flameholder. 

The  vortices  grow  as  they  move  downstream  and,  after  a  fluid  dynamic 
time  delay  If,  their  interaction  with  the  lower  wall  produces  an 
increased  mixing  between  the  unburnt  mixture  and  the  completely 
burnt  gas  in  the  recirculations  zone.  After  a  further  time  delay  or 
chemical  induction  time  Ic,  which  is  fixed  primarily  by  chemcial 
parameters,  this  increase  in  mixing  rate  will  produce  an  increase  in 
the  heat  release  rate. 


When  the  total  delay  time,  I  f  +  Ic,  has  the  appropriate 
magnitude  with  respect  to  the  period  of  the  oscillation,  this 
■fluctuation  in  the  heat  release  will  dr  i  ve  the  pressure  oscillation. 
If,  however,  the  total  delay  is  too  long,  the  fluctuation  will 
extract  energy  from  the  acoustic  field  and  the  pressure  oscillations 
will  decay . 

This  picture  is  over  simplified  primarily  becasue  the  heat 
addition  fluctuation  occurs  over  a  period  which  is  comparable  to 
that  of  the  oscillation  and  is  spread  out  over  the  whole  combustion 
chamber.  Hence,  to  discern  the  overall  effect  of  the  oscillation, 
it  is  necessary  to  use  a  representation  of  the  Rayleigh  criterion 
which  involves  an  integral  over  a  whole  period  of  the  oscillation 
and  which  extends  over  the  entire  volume  in  which  heat  is  added  to 
the  gas  stream. 

Frequency  Estimates-  The  acoustic  modes  are  determined  by 
the  longitudinal  geometry  of  the  entire  system  including  the  supply 
chambers  as  well  as  the  combustor.  A  one-dimensional  acoustic 
analysis  of  the  system  has  been  used  to  determine  the  natural 
frequencies  and  the  corresponding  mode  shapes.  These  have  been 
compared  to  the  experimental  frequencies  of  oscillation  and  the 
experimental  pressure  mode  shapes  when  operating  conditions  yield 
large  pressure  oscillations.  The  frequencies  agree  quite  well  with 
the  experiment  and  the  theory  predicts  the  location  of  modal  nodes 
and  the  mode  shapes  fairly  accurately.  The  frequency  of  oscillation 
depends  in  a  predictable  manner  upon  the  geometry  of  the 
flameholder,  the  mean  flow  speed,  the  fuel  type,  and  the  equivalence 
ratio. 


Influence  of  Chemical  Parameters-  During  the  past  year,  a 
parametric  study  was  performed  to  investigate  the  dependence  of  the 
vortex  shedding  frequency  upon  variations  in  the  mean  flow  speed  and 
equivalence  ratio  for  a  fixed  geometry  and  fuel  type.  The  results 
reveal  that  small  changes  in  equivalence  ratio  can  lead  to 
excitation  of  acoustic  mooes  of  greatly  varying  frequency.  Thus,  the 
chemical  induction  time  was  determined  to  be  important  even  in  low 
frequency  instabilities. 


The  chemical  time  dependence  on  the  parameters  of  the  system  is 
consistent  with  the  model  described  above.  For  example,  the  model 
indicates  that  shorter  time  delays  allow  the  mechanism  to  support 
higher  frequency  oscillations.  The  chemical  induction  time  is 
decreased  as  the  fuel-air  ratio  approaches  the  stoichiometric  value 
or  when  the  fraction  of  hydrogen  in  the  me  thane-hydrogen  fuel 
mixture  is  increased.  In  agreement  with  these  expectations,  we 
observe  in  the  experiments  that  high  frequency  modes  can  be  excited 
when  we  use  fuel-air  ratios  near  stoichiometric  but  not  for  lean  or 
rich  mixtures.  Similarly,  when  fuels  which  contain  high  fractions 
of  hydrogen  are  used,  high  frequency  modes  are  excited. 
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Fluid  Dynamic  Parameters-  Reduction  of  the  step  height 
results  in  the  excitation  of  higher  -frequency  modes  as  we  expect 
because  the  -fluid  dynamic  delay  time  is  reduced.  In  addition  this 
result  supports  the  view  that  the  influence  of  mixing  as  the  vortex 
structure  impinges  on  the  lower  wall  of  the  combustor  is  important 
for  most  instabilities  observed. 

We  have  also  found  experimentally  that  pressure  amplitudes 
increase  when  the  frequency  of  the  osciallations  increases;  this 
result  is  in  agreement  with  the  cal cualt ions  carried  out  in  this 
program  by  Hendricks  who  showed  that  the  fluid  dynamic  time  delay 
decreased  almost  lineraly  with  the  amplitude  of  the  pressure 
fluctuation. 

Finally,  higher  frequency  modes  are  observed  at  higher  mean 
r i ow  speeds.  This  effect  is  not  well  understood  but  may  be  caused 
by  decreases  in  the  mixture  fraction  of  cold  reactants  at  high 
speeds  resulting  in  a  shorter  chemical  reaction  time. 

Stability  Criteria-  One  of  the  primary  results  of  the 
research  carried  out  during  the  past  year  was  the  collection  of 
complete  sets  of  light  intensity  and  pressure  fluctuations  data  for 
the  whole  duct  and  for  a  number  of  different  configurations  of  the 
system.  The  light  intensity  data  are  used  to  determine  the  heat 
addition  fluctuations,  and  these,  in  conjunction  with  pressure 
fluctuation  data,  have  been  used  to  determine  the  input  of  energy  to 
an  acoustic  mode  by  the  heat  addition  process. 

Comparison  of  frames  from  shadowgraph  movies  with  measured 
values  of  the  instantaneous  distributions  of  the  radiation  intensity 
in  the  combustion  chamber  during  combustion  instability  reveal  that 
the  volumetric  rate  of  heat  release  due  to  combustion  in  the  regions 
near  the  flameholder  continue  to  rise  long  after  the  large  vortex 
structure  passes.  At  locations  farther  down  stream,  the  amplitude 
of  the  fluctuating  radiation  intensity  signal  is  larger  and  the  peak 
in  reaction  occurs  only  shortly  after  the  vortex  passes.  Thus,  an 
apparent  combustion  wave  is  produced  in  the  duct  by  the  passage  of 
the  vortex  and  the  front  of  this  wave  moves  at  a  nearly  constant 
velocity  that  is  about  twice  the  mean  flow  speed  for  standard 
operating  conditions.  The  accelerating  vortex  structure,  on  the 
other  hand,  was  found  to  move  at  about  half  the  mean  flow  speed 
measured  at  the  flameholder  lip. 

Data  for  the  distribution  of  combustion  intensity  versus 
distance  downstream  form  the  flame  holder  are  shown  in  Fig.  2,  3 
for  eight  different  times  during  a  complete  cycle;  correspond i ng 
shadowgraph  photographs  for  same  time  intervals  are  also  shown. 

The  combustion  wave  described  above  is  clearly  visible. 


3.  Vortex  Shedding  at  188  hz ;  Second  Halt-Cycle 

The  role  of  Rayleigh's  Cr , t ?r i on  was  examined  in  detail  as 
the  predominant  process  driving  and  damping  instabilities  of 
combustion  chambers.  Energy  considerations  show  that  the  potential 
for  augmentating  the  acoustic  field  by  coupling  of  pressure 
fluctuations  and  heat  release  rate  fluctuations  is  about  two  orders 
of  magnitude  greater  than  other  energy  transfer  mechanisms.  This 
observation  suggested  measuring  phase  difference  between  the 
fluctuating  pressure  and  the  oscillatory  radiation  intensity  signal 
at  locations  throughout  the  combustion  chamber  for  several  example 
operating  conditions. 


The  results  indicate  that  regions  of  strong  damping  general  )v 
occur  immediately  upstream  of  the  location  at  which  the  amplitude  of 
the  radiation  intensity  is  a  peak.  Similarly,  regions  ot  strong 
damping  ot  the  acoustic  energy  usually  occur  just  downstream  ot  ♦’his 
location.  Thus,  the  location  ot  the  largest  heat  release  occurs  at 
a  time  during  the  pressure  oscillation  that  results  in  no  driving  or 
damp i ng  ot  the  acoust i c  t i e 1 d  after  a  finite  limit  cycle  has  been 
reached . 

Modelling  the  feedback  me  chan  ism  to  determine  the  phase 
difference  between  the  pressure  and  the  heat  release  oscillations 
traditionally  uses  a  pressure-sensitive  time  delay  to  predict  linear 
stability  char ac ter i st i cs .  This  has  given  considerable  success  for 
liquid-propellant  and  solid-propellant  systems.  However,  for 
combustion  systems  in  which  the  flame  is  stabilized  in  a 
recirculation  zone,  the  feedback  mechanism  associated  with  vortex 
shedding  phenomena  is  sensitive  to  the  velocity  fluctuations  at  the 
f 1 ameho 1 der . 

A  linear  stability  mode  1  based  on  a  velocity-sensitive  time  lag 
is  used  for  the  geometry  of  the  laboratory  combustor.  The 
combustion  is  model  led  by  a  volumetric  source  downstream  of  the 
flameholder  whose  amplitude  is  related  to  the  amplitude  of  the 
velocity  fluctuations  over  the  flameholder  and  whose  phase  lags  the 
same  velocity  by  some  time  delay.  The  solution  method  requires  that 
particular  values  of  the  time  delay  and  the  interaction  index 
determine  the  amplification  rate  and  frequency  of  the  oscillation. 
For  feedback  within  one  cycle  of  oscillation  the  time  delay  required 
for  instability  is  generally  shorter  for  higher  frequencies.  Large 
changes  in  the  phase  difference  between  the  pressure  fluctuation  and 
the  velocity  fluctuations  occur  near  natural  frequencies  of  the 
system.  This  has  a  strong  influence  on  the  stability  boundaries 
predicted  for  the  system. 

Finally,  possible  mechanisms  of  nonlinear  growth  are  described 
within  the  framework  of  this  model.  Experimental  results  for 
particular  operating  conditions  reveal  that  a  combination  of 
mechanisms  that  are  discussed  may  be  responsible  for  limit  cycle 
behavior.  Ex per i men  t  a  1  results  show  that  the  amp litudes  of  high 
frequency  modes  are  generally  larger  than  the  amplitudes  of  low 
frequency  modes  and  that  the  frequency  of  oscillation  is  usually 
higher  than  predicted.  These  results  suggest  a 

nonlinear  evolution  that  involves  an  increase  in  the  frequency  of 
oscillation  and  a  decrease  in  the  time  delay  parameter  as  the 
amp litude  of  the  oscillation  reaches  a  finite-limiting  amp  1  i  t ude . 


The  primary  work  carried  out  under  this  grant  has  been  the 
design  and  construction  of  a  new  test  duct  tor  the  proposed 
experiments,  uur  aim  was  to  develop  an  apparatus  in  which  we  could 
produce  vortices  at  the  inter-face  between  burnt  gas  and  a  fresh  fuel 
air  mixture  and  to  observe  the  influence  of  various  fluid  dynamic 
and  chemical  parameters  on  the  rate  of  combustion  in  these  vortices. 

One  of  the  parameters  we  want  to  observe  is  the  light  emitted 
during  the  combustion  process  which  is  a  good  measure  of  the  local 
volumetric  heat  release  by  the  combustion  process.  Ule  shall  do  this 
with  a  camera  based  on  a  charge  injected  device  which  has  an  array 
of  380  by  240  pixels  and  employs  a  gated  image  intensification 
system  which  is  requried  to  sense  the  low  light  levels  produced  by 
the  combustion  of  methane  and  air.  This  system  will  allow  us  to 
observe  a  region  of  about  15  cm  by  24  cm  with  useful  resolution  it 
requires  about  100,000  data  points  per  picture.  The  data  acquisition 
rate  required  to  follow  the  deve lop erne nt  of  a  vortex  with  this 
camera  is  over  100  million  per  second  and  is  larger  than  our  system 
can  accomodate.  Hence,  we  shall  observe  a  periodic  sequence  of 
vortices  and  use  a  phase  averaging  technique. 

Based  on  our  experience  in  the  combustion  instability 
experiment,  we  have  chosen  to  use  the  same  rearward  facing  step 
geometry  to  produce  the  vortices.  However,  we  have  modified  the 
combustion  duct  so  that  we  may  control  the  frequency  and  growth  rate 
of  the  vortices.  In  addition  we  have  increased  the  depth  of  the  test 
section  so  that  the  vortices  can  grow  to  a  much  1 arger  scale  before 
being  influenced  by  the  lower  wall  of  the  combustion  chamber. 

New  Test  Facility-  A  schematic  diagram  of  the  new  test 
facility  is  shown  in  Fig. 4.  The  most  evident  change  in  the 
combustion  chamber  is  the  use  of  a  duct  with  a  height  that  can  be 
varied  bewtween  2.5  and  15  cm.  The  maximum  depth  will  allow  vortices 
to  reach  scales  of  10  cm  before  the  inteference  effects  of  the  lower 
wall  becomes  important. 

In  the  previous  experiments,  the  frequency  of  the  oscillation 
selected  by  the  instability  was  one  of  the  acoustic  modes  of  the 
system  and  the  mode  selected  was  a  complex  function  of  the  length  of 
the  plenuim  chamber,  the  geometry  of  the  rest  of  the  system  and  the 
chemical  parameters  of  the  fuel-air  mixture.  In  the  present 
apparatus,  we  will  use  this  result  to  our  advantage  by  controlling 
the  length  of  the  plenum  chamber  to  fix  the  frequency  of  one  of  the 
acoustic  modes  to  match  the  desired  frequency.  We  gain  further 
control  over  the  frequency  of  the  vortex  shedding  by  using  a  siren, 
placed  at  the  upstream  end  of  the  combustion  chamber,  to  excite  the 
duct  at  the  desired  frequency.  The  amplitude  of  the  energy  supplied 
by  the  siren  to  the  desired  mode  can  be  varied  by  passing  part  of 
the  fuel -air  mixture  through  the  siren  and  part  through  a  bypass 
around  the  siren. 
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Diagram  of  New  Experimental  Facility 


This  system  will  allow  us  to  control  the  amplitude  and 
frequency  of  the  velocity  and  pressure  fluctuations  at  the  flame 
holder  lip.  Thus,  we  will  be  able  to  control  the  frequency  of  vortex 
shedding  from  the  lip.  Futher,  our  previous  work  has  shown  that  the 
rate  of  growth  of  the  vortices  is  a  function  of  the  amplitude  of  the 
fluctuations  at  the  flame  holder  lip.  The  bypass  arrangement  which 
we  have  provided  for  mixture  injection  allows  us  active  control  of 
the  amplitude  as  well  as  the  frequency  of  the  fluctuations  and 
thereby  control  the  rate  of  growth  of  the  vortices.  The  apparatus 
will  produce  vortices  at  frequencies  between  100  to  600  Hz  and  will 
allow  us  to  observe  vortices  with  scales  up  to  15  cm.  Construction 
has  been  completed  and  initial  tests  are  in  preparation. 


I  ns t rumen t a  t i on-  The  primary  new  item  of  instrumentation  is 
the  gated  image  intensified  CID  camera  which  will  be  used  to  measure 
the  light  emitted  by  the  combustion  process  and  thus  to  estimate  the 
rate  of  heat  release  throughout  the  vortex.  Photographs  will  be 
obtained  at  a  frequency  of  30  frames  per  second  and  will  be  stored 
an  a  MCR.  Later,  individual  frames  will  be  selected  from  this  video 
movie  with  a  frame  grabber  and  the  intensity  data  will  be  digitized 
for  further  analysis. 
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In  addi t i on ,  we  wi 1 1  use  a  new  LDU  system  to  measure  one  op  two 
velocity  components  throughout  the  flow  field.  This  Disa  system  with 
a  single  velocity  component  capability  was  acquired  during  the 
1986-67  contract  year  and  the  two  component  system  will  be  completed 
early  in  the  1987-88  contract  year. 

Other  parameters  measured  will  include:  the  velocity 
■fluctutation  at  the  flame  holder  lip,  obtained  with  a  conventional 
hot  wire  system;  pressure  and  pressure  fluctuations  on  one  wall  of 
the  duct;  single  frame  shadowgraph  pictures;  and  8,000  frames  per 
second  shadowgraph  movies.  Finally,  we  have  developed  an  ion  probe 
which  will  allow  us  to  determine  the  ion  density  in  a  volume  of  gas 
corresponding  to  a  1  mil  lime  ter  cube.  Because  the  regions  where 
combustion  is  taking  place  are  also  regions  in  which  the  gas  is 
highly  ionized,  this  i nsturmen  t  will  be  used  in  conjunction  with  the 
flow  visualization  techniques  to  make  certain  that  we  can  define  the 
regions  of  active  combustion  accurately. 


7.  Two  Inter-facial  Contours 
Developed  Within  Annulus 


6.  Strain  and  Rotation  o-f 
Element  Initially  on 
Horizontal  Axis 


between  two  concentric  circles,  Fig. 6.  The  element  makes  an  angle 
j {  with  respect  to  the  local  tangent  to  the  circle  and 

~  57  •  3  A  3 

By  i  n  tegrat  i  ng  ,  from  Eq  .  2,  we  obtain  the  remarkably  simple 
result 


A  **•  Vjt/r  4 

The  length  o-f  the  element  o-f  interface  is 

<u*  --  s 

The  motion  reduces  with  time,  the  interface  is  continuously 

being  stretched  along  its  length,  and  at  some  time  later  makes 
a  complete  loop  about  the  origin,  Fig.  7.  When  these  spirals  are 
closely  spaced,  this  angle  is  nearly  constant  and  Eq .  5  may  be 
integrated  about  a  complete  circle  giving 

zvr  =.  (/  +  6 

where  -  r,  is  the  radial  extent  of  the  circular  annulus  in  which 
the  angle  of  the  deformed  interface  is  ^  .  In  actuality  this 

strip  contains  two  interfacial  elements,  Fig.  4,  one  originating  on 
the  right  horizontal  axis,  the  other  on  the  left.  One  may  interpret 


-13- 


Eq .  6  as  giving  the  width  of  an  annulus  containing  a  single 
stretched  line  element  circling  the  origin, 


7Tr  /  (  /  +  al) 


It  is  convenient  to  construct  a  local  coordinate  system,  Fig. 8,  and 
to  express  the  local  -flow  -field  in  these  coordinates.  The  local  x 
and  y  velocity  components  are 


u  -  I  (i.  *~~i.*0*  +  ( 

v~  -  -  (-z  ?A/di 


The  inter-face  is  always  alligned  with  the  moving  x 
hence  the  straining  rate  along  the  interface  is 


-  axis  and 


£  =  (i:  r  A/(i4A1)% 


8  Coordinate  System  for  Local  Analysis 


Fast  Chemistry  in  a  Cortex  Field-  By  "fast  chemical 
reaction  rates"  it  is  implied  that  the  two  reactants  cannot  coexist 
on  the  time  scale  of  interest  in  the  problem.  Physically,  this 
requires  that  the  chemical  time  be  short  in  comparison  to  the  time 
associated  with  diffusive  mixing.  This  -feature  introduces  two 
significant  features  into  the  analysis:  i)  There  is  neglegible 
effect  of  remote  strips  of  reactant  upon  the  one  under- 
consideration,  and  ii)  the  process  is  complete  in  a  given  strip  when 
the  reactants  are  consumed,  which  occurs  in  a  finite  time. 

For  the  local  analysis,  therefore,  it  is  sufficient  to  consider 
a  single  in  an  infinite  plane  domain  with  k,  initially  in  the  upper 
half-plane  and  in  the  lower  half-plane.  As  time  progresses, 

the  two  reactants  diffuse  toward  their  interface,  at  which  K,  and 
both  vanish  because  of  fast  chemistry,  and  react  in  their 
stoichiometric  ratio  to  produce  the  product.  The  product,  in  turn, 
diffuses  back  into  the  upper  and  lower  half— planes.  If  there  were  no 
strain  rate,  the  reactant  distribution  would  be 

K(j  ±  erf  I  7~/( 

where  a  sto i ch i ome tr i c  ratio  of  unity  has  been  assumed.  The  rate  at 
which  the  reactant  Ki  is  consumed  at  the  reaction  interface  is 

-  (»/vt)yx  12 

On  the  other  hand,  when  the  local  interfacial  region  is  being 
strained,  the  diffusion  equation  takes  a  more  complex  form. 

The  appropriate  solution  for  K,  is  then 

K ,  --  erf[i/(^T)^j  13 

and  the  consumption  rate  of  reactant  1 

iS  i 

r 

v  7T  r  ; 

Recalling  Eq.  10  giving  the  strain  rate  £(4)  in  the  vortex, 

the  rate  at  which  the  reactant  is  being  consumed  per  unit  area 

of  i n  ter f ace  is 


per  unit  area  of  interface 
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Now  an  element  of  interface  originally  of  length  dr  has  been 
stretched  into  a  length  d-  S  ,  Eq .  5,  an  hence  the  reactant 
consumption  rate  by  this  element  is 

r  ft  y'/*- 

V(r,t)JL4  .  C%) 4  |  j  + 
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In  the  unreactec  -  '.cte  an  annulus  ATiret-r  contains  halt  reactant  1 
with  concentration  «,  -  <  ,  and  hal-f  by  reactant  2  with 

concentration  .  Moreover,  as  was  pointed  out  in  connection 

with  Fig.  7,  there  are  two  inter-faces  within  this  annulus.  Thus 
the  reactant  consumption  rate,  Eq .  1*,  integrated  over  time,  will 
at  some  time  t*  consume  all  of  the  reactant  1  within  the  annulus 
Expressed  anal y t i cal  I y , 

Z  f  (Vt'i  idsjut  -  7?  rjtr 
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This  relation  may  be  written  in  detail  using  Eq .  16.  The  resulting 
expression  may  be  cast  in  the  -form 

s 


where 


hit  - 


P  /l  v  /s  /n  t  * 


Because  ft#  i  s  a  number  ^  ^  ,  Eq  .  18  gives  the  reacted  core 

radius  at  any  time  t  after  the  start  o-f  motion. 


The  result  is  independent  ot  Reynolds  number  except  tor  very 
large  values  of  yK-*  ■  The  asymptotic  behavior  tor  large  is 

f  n  / (  P l/i  ov*  t  )  -  o.  S' o*f 

22 

The  corresponding  result  tor  low  values  ot  is 


^  /( r V<  V7*2-  -•  0.0  23 

This  result  tor  low  values  ot  is  new  and  of  importance  under 
circumstances  where  the  viscous  dittusion  and  the  dittusion  due  to 
distortion  are  ot  the  same  order  ot  magnitude.  In  tact  the  new 
dimensionless  ratio  ,  given  in  Eq.21,  represents  the  ratio  ot 

the  rate  ot  growth  ot  the  burned  core  radius  to  the  rate  ot  growth 
ot  the  viscous  core. 
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No  t  a  t i on  - 

D  Coefficient  of  molecular  dit' -fusion 

dr  Interfacial  element  length  before  distortion 

ci S  Interfacial  element  length  after  distortion 

I  Defined  in  Eq.  1? 

A  &  Polar  coordinates  for  vortex  flow 
A1*  Radius  of  burned  core 

R.  Reynolds  number  of  vortex, 

i  Time 

£ #  Defined  by  Eq .  17 

W ,  V/^  Radial  and  tangential  velocity  components 
V  Reactant  consumption  rate  per  unit  area  of  interface 

i  +  A 


r/zTT? 


1 

7 


Local  cartesian  coordinates  of  interfacial  element 
j  C  i  A  )dt 


P> 

r 

e 


Defined  in  Eq .  20 
Circulation  of  vortex 
Element  strain  rate,  Eq .  10 
X-,  t  Concentration  of  Reactant  1  and  2  respectively 
A  Def i ned  by  Eq .  4 

yUe »,  Ratio  of  diffusive  spreading  to  viscous  spreading, 

Eq.  21 

9  Kinematic  viscosity 

¥  Angle  of  interfacial  element,  Eq .  3 
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